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Abstract 



A recent experiment by Ma and collaborators shows that the gap 
anisotropy in Bi2Sr2CaiCu208+x is strongly temperature dependent. In par- 
ticular, the superconducting gap along the T — M direction shows a weaker 
temperature dependence than the gap along the T — X direction which de- 
creases rapidly with temperature. We explain this novel feature as a natural 
consequence of the interlayer tunneling mechanism of superconductivity. 
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The nature of the order parameter in the cuprate superconductors is of central importance 
to the study of high temperature superconductivity. While several experiments such as 
IR reflectivity ||1||, Raman spectroscopy 0, NMR and Josephson tunneling |Q imply 
that a conventional s-wave order parameter as in BCS theory is unlikely, Angle Resolved 
Photoemission Spectroscopy (ARPES) provides direct evidence that the gap in the cuprate 
superconductors is highly anisotropic The main results that can be inferred from low 
temperature ARPES |^ are : (i) the superconducting gap attains its maximal value along 
the r — M direction (Cu-0 bond direction in real space); (ii) the gap is smallest along the 
r — X{Y) direction (diagonal to the Cu-0 bond direction) and (iii) there is a monotonic 
increase in the magnitude of the gap from its smallest value along the T — X line to the 
largest value along F — M, thereby indicating true anisotropy. 

In a recent paper 0, Ma and co-workers have presented the first detailed analysis of 
photoemission spectra of superconducting Bi 2212 observed along T—M and F— X directions 
at different temperatures. Their results show that the gap anisotropy observed in this 
material is strongly dependent on temperature, contrary to what happens in conventional 
anisotropic superconductors such as Pb. In particular they show that the gaps observed 
along the two high symmetry directions obey different temperature dependences. The gap 
along the T — M direction shows a very weak temperature dependence whereas the gap 
along T — X decreases rapidly as temperature increases. Consequently, the gap anisotropy 
increases with temperature by about a factor of 8 before falling to zero at Tc . Since these 
results reflect the property of the superconducting condensate in the two different high 
symmetry directions, they impose constraints on possible mechanisms of high temperature 
superconductivity. 

In this letter, we show that these results can be explained by the interlayer tunneling 
mechanism of high temperature superconductivity . We show that the gap equation result- 
ing from interlayer tunneling ^ readily distinguishes between the gaps along the T — M and 
F — X directions and the temperature dependences of the gaps along these two high symme- 
try directions are completely different. Therefore, a temperature dependent gap anisotropy 

2 



follows very naturally from the interlayer tunneling mechanism. Our analysis shows that 
the gap anisotropy increases rapidly with temperature and the results are consistent with 
those reported in ref. 

We begin by writing the gap equation from interlayer tunneling 0, 

= '""^^ + ^' '""^^ • 

This equation can be obtained by considering two close Cu-0 layers as in Bi 2212 coupled 
by a Josephson tunneling term of the form, 

k 

where t is a band structure parameter in the dispersion of electrons along the Cu-0 plane 
and tj_{k) is the bare single electron hopping term between the two coupled layers c and d. 
The quantity Tj{k) in the right hand side of equation (1) is given by Tj{k) = iti^. The 
dispersion of electrons along the Cu-0 plane is chosen to be of the form 

e{k) = — 2t {coskx + cosky) + At' cosk^ cosky , 

with t = 0.25 eV and j = 0.45. We also choose ep = -0.45 eV corresponding to a Fermi 
surface which is closed around the F - point. These choices are inspired by band structure 



calculations |T0|. Note that the Josephson coupling term in Hj conserves the individual 
momenta of the electrons that get paired by hopping across the coupled layers. This is as 
opposed to a BCS scattering term which would only conserve the center of mass momenta 
of the pairs. This is the origin of all features that are unique to the interlayer tunneling 
mechanism. The second term in the right hand side of equation (1) is obtained by postulating 
that the dominant in-plane pairing mechanism is the electron-phonon interaction. The 
primed sum in equation (1) is over a shell about the Fermi surface of width huD, the Debye 
energy. 

Let us first consider equation (1) in the limit of zero temperature 0. We then have 

Afc = + , (2) 
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where A'' = Vbcs J2'q is a finite s-wave component of tlie gap. Obviously, tlie anisotropy 
in tlie gap at zero temperatures is principally due to the momentum dependence of Tj{k), 
viz., that of t±{k). This momentum dependence can be inferred from electronic structure 
calculations. As shown in ref. it is adequate to choose Tj{k) = ^(cosA;^; — cosky)'^ to 
reproduce the results of such calculations. With this choice then, it is obvious from equation 
(2) that the gap is smallest when Tj{k) = (as long as Tj > Vbcs)- This happens when 
kx = ky. Therefore along the T — X{Y) directions where k^ = ky, the gap attains its smallest 
value A*. The largest value of the gap is obtained when the Fermi surface includes the points 
where Tj{k) is maximum, viz.,(0, ±7r) and (±7r, 0). In this case, the value of the gap is given 
by ^ + A'' as can be seen from equation (2). So, we see that the zero temperature gap 
from interlayer tunneling is highly anisotropic with maxima at (0, ±7r) and (±7r, 0) and 
minima along the k^ = ky line. These results are in agreement with the ARPES data on 
gap anisotropy at low temperatures. We now show that this gap anisotropy increases as 
temperature increases. 

To see this we go back to the gap equation at finite temperature, equation (1). We first 
consider the gap at (0, vr) on the T — M line, A(r — M) which is given by 

A(r - M) = ^ tanh ^^^^^~ + A%T) , (3) 



where 



A%T) = Vbcs E'TT^tanh 
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The temperature dependence of A(r — M) will be governed by the first term in the right hand 
side of equation (3) since we always choose to work in the limit of A^{T) << Tj. This limit 
corresponds to the physical choice of Tj giving rise to high transition temperatures rather 
than the electron-phonon interaction Vbcs- Since the structure of the gap equation (3) is 
very different from that of the BCS gap equation, it is clear that the temperature dependence 
of A(r — M) will be unlike that of a BCS gap. For instance, A(r — M) will fall steeper near 



Tc than a BCS gap would. This is seen most simply in the limit A^{T) ||rT|. In this 
limit, the zero temperature gap Ao(r — M) = It is readily seen that the temperature at 



T 

which A(r — M) falls to half its zero-temperature value is given by — ^'^^ typical 
values of t±_ and t, this temperature corresponds to T ~ 0.95Tc. Consequently, we expect 
A(r — M) to show a weak temperature dependence at low and intermediate temperatures. 
It should be emphasized that this temperature dependence (as given by the first term in 
the right hand side of equation (1)) is directly due to the "momentum-space locality" of the 
Josephson interaction Hj. 

We have solved equation (3) for temperatures ranging from to T^. Our choice of 
parameters are : t = 0.25 eV, t' = 0.1125 eV and ep = -0.45 ev as mentioned earlier, t± = 
0.091 eV, Vbcs = 0.06 eV and hiOc = 0.02 eV. This choice of parameters leads to a purely 
in-plane of ~ 5K and a bulk of 83K. The zero temperature gap at (0, tt), Ao(r — M) is 
found to be 18.09 meV. In fig.(l), we have shown our results for the temperature dependence 
of A(r — M) and compared them with the experimental results of Ma et al 0. Note in 
particular that the measured gap (a) is very weakly temperature dependent at low and 
intermediate temperatures and (b) falls very steeply near T^. As we mentioned earlier, both 
these features follow directly from equation (3) as consequences of the Josephson interaction 
Hj. In view of this, we suggest that PES along the T — M direction actually probes the gap 
resulting from interlayer tunneling. 

We now consider the gap along the F — X{Y) directions where = ky. In this case, the 
first term in the right hand side of equation (1) drops out and the gap A(r — X) is given by 

AiT-X) = A^{T) = FBc^E'^tanh^ . (4) 

It is obvious that the temperature dependence of A(r — X) will be different from that of 
A(r — M). Note that equation (4) looks very much like the BCS gap equation. However 
it should be emphasized that A(r — X) is not a BCS gap since the sum in equation (4) 
also contains the enhancement effects of the Josephson interaction. In fact our results 
show that the gap A(r — X) falls much faster from its zero temperature value than a 
BCS gap does. Therefore we get a non BCS-like temperature dependence for A(r — X) 
as well. It is not necessary to solve equation (4) to obtain the temperature dependence of 



A(r — X) since A(r — X) = and we have already solved for A*(r) when we obtained 
the temperature dependence of A(r — M). From our results for A(r — M), it is easy to see 
why A(r — X) decreases rapidly with temperature. For instance, as T increases from O.lTc 
to 0.8Tc, A(r — M) decreases from 18 meV to 13 meV, we find that A* has to decrease by 
a factor of 3 for equation (3) to be self consistent. This is because of the first term in the 
right hand side of equation (3) which depends weakly on temperature. Therefore, the rapid 
decrease of A(r — X) with temperature is actually a consequence of the weak temperature 
dependence of A(r— M). We have obtained A(r— X) for various temperatures with the same 
set of parameters as before. We find that there is a quantitative discrepancy between our 
results and those of Ma et al. With our choice of parameters, we find the zero temperature 
value of the gap along F — X{Y), Ao(F — X) ~ 2 meV. This value is closer to the values 
reported in earlier ARPES experiments 0. On the other hand, the low temperature value 
quoted by Ma et al is A(F - X) = 10±2 meV at T = 0.48 T^. While it is possible for us to 
tune Vbcs aiid obtain larger values of A(F — X) without altering Tj and substantially, 
we have not done so for the following reasons. The first is that the observed values of the 
gaps along F — X are very sensitive to sample quality and the time elapsed between cleaving 
and observation of the spectrum Secondly, we feel that given the level of approximation 
involved in this modelling, such a fine-tuning of parameters is unwarranted. Instead, to 
see how the gap anisotropy increases with temperature, we only consider the temperature 
dependence of A(F — X) normalized to its zero temperature value. This is still meaningful 
as we are only interested in seeing how the gap anisotropy grows from its low temperature 
value and not in the absolute values of the gaps themselves. On comparing our results for 
the temperature dependence of A(F — X), with those of Ma et al, we find that at low and 
intermediate temperatures there is a discrepancy of 5-10% . At temperatures close to T^, 
the discrepancy is slightly more. This is because the experimental results indicate that the 
gap becomes vanishingly small at 0.84 whereas a mean field theory such as ours would 
always produce gaps A(F — M) and A(F — X) that vanish self consistently at identical 
temperatures. 
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In fig. (2), we have shown the temperature dependences of both A(r — M) and A(r — X) 
as obtained from equation (3). The solid hne shows the temperature dependence of the 
former and the dashed line shows that of the latter. We find that A(r — X) falls to half its 
zero-temperature value at temperatures as low as O.QTJ The observed experimental value 
of this temperature is 0.64Tc Recall that a conventional BCS gap falls to half its zero- 
temperature value at T ~ O.OTc. As T — Tc, we find that A(r — X) decreases by an order 
of magnitude at T ~ 0.9Tc while at this temperature, A(r — M) has only decreased by a 
factor of 2. Thus, the gap anisotropy increases from its low temperature value by a factor 
of 5. 

For completeness, we have also investigated the temperature dependence of the gap 
anisotropy when the order parameter has mixed s- and d- wave symmetries. We have done 
this because the low temperature ARPES results can also be explained by any model 
exhibiting d-wave superconductivity However the results of ref. 0, particularly the 

temperature dependent gap anisotropy cannot be obtained from an order parameter having 
pure d-wave symmetry. This is because the anisotropy ratio of the gap for any two k 
vectors is temperature independent if the gap function has a purely d-wave character. One 
possibility is an order parameter with mixed s- and d- symmetries. The pairing potential 
leading to such an order parameter will be of the form 

Vkk' = Vq + V^i(cosA;2,. — cosky) (cosA;^,.' — coskyi) . 

We have solved the gap equation resulting from this interaction for various choices of Vq and 
Vi without including the interlayer Josephson interaction. We find that at low temperatures, 
the gap anisotropy always decreases irrespective of the relative strengths of Vq and Vi. At 
intermediate temperatures, it is possible to obtain an increase in the gap anisotropy by 
fine-tuning the parameters but this increase is only marginal. These results underscore 
the importance of the Josephson interaction Hj in establishing and enhancing the gap 
anisotropy. 

Finally, we address the question of the sensitivity of our results vis a vis our choice of 
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parameters. We have solved equation (1) for various choices of t± and Vbcs keeping the 
other parameters fixed. We find that as long as the value of the s-wave component of the 
gap is much smaller than Tj, i.e., in the limit of interlayer tunneling being stronger than any 
intralayer interaction, the gap anisotropy always increases with temperature and differences 
in results are only quantitative. 

To conclude, we have shown that the interlayer tunneling mechanism produces a gap 
anisotropy that grows with temperature. This is because the gap equation from interlayer 
tunneling leads to different temperature dependences for the gaps along the two high sym- 
metry directions in Bi 2212. The gap along the T — M direction shows a much weaker 
temperature dependence than the gap along the T — X direction which decreases rapidly as 
temperature increases. Consequently, the gap anisotropy increases with temperature by a 
factor of ~ 5. This is in good agreement with experimental results that show an increase 
by a factor of ~ 8. Our results and those of Chakravarty et al. show that the interlayer 
tunneling mechanism of high temperature superconductivity can account satisfactorily for 
the ARPES experiments in superconducting Bi 2212. 
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Figure Captions 

All the results were obtained with the following choice of parameters : t = 0.25 eV, t' = 
0.1125 eV, ep = -0.45 eV, t± = 0.091 eV, Vbcs = 0.06 eV and Hud = 0.02 eV. 

1. A(r — M) as a function of the reduced temperature ^. Solid line is as calculated from 
equation (3) in text. The experimental data is from ref. |^. 

2. Temperature dependence of the gaps (normalized to zero temperature values) along 
r — M (solid hne) and T — X (dashed line) as a function of reduced temperature 
obtained from equation (3). 
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